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Abstract 

We consider conductivities of two-dimensional lattice electrons in a magnetic 
field. We focus on systems where the flux per plaquette (ft is irrational (in- 
commensurate flux). To realize the system with the incommensurate flux, we 
consider a series of systems with commensurate fluxes which converge to the 
irrational value. We have calculated a real part of the longitudinal conductiv- 
ity a xx (uj). Using a scaling analysis, we have found ^a xx {ijj) behaves as 1/w 7 
(7 = 0.55) when (ft = r, (r = y/ ^~ 1 ) and the Fermi energy is near zero. This 
behavior is closely related to the known scaling behavior of the spectrum. 
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I. INTRODUCTION 



The electronic properties of the two-dimensional periodic systems in a magnetic field have 
been studied extensively. Even for non-interacting electrons, the various physical quantities 
(for example, the wavefunctions, and the energy spectra) exhibit extremely rich behaviors 
P|-|5|,|S| jr5| , P^7[] and it has been attracted great attentions in relation to the quantum Hall effect 
|7HT2], one-dimensional quasiperiodic systems ]T3— [T7[| , flux states for the high-T c supercon- 
ductivity [P^-f2T|. The algebraic structure of this problem has also been revealed recently 



Consider the tight-binding Hamiltonian on the square lattice 

H° = - e^cjcj + h.c, (1.1) 

<ij> 

where cj(cj) is a creation (annihilation) operator of an electron at site i. The summation 
is taken over the nearest neighbor sites. The hopping amplitude between link (ij) is set to 
be unity. The phase factor 0^ is defined on link (ij). The magnetic flux per plaquette is 
(ft = ^ I]pi aque tte &ij i n units of magnetic flux quantum he/ e. 

When <p is rational, i.e. = p/q with mutual prime integers p and q, the spectrum 
consists of q bands with finite widths. The wave functions are the extended Bloch functions. 
Many interesting phenomena related to the quantum Hall effects are discussed |8Hl2|. 



When <p is irrational (incommensurate), the system exhibits novel structures ||. It is 
known that the electronic state and the energy spectrum have various singular natures. The 
energy spectrum is a Cartor set which consists of infinitely many bands with zero width 
T3| , |T5|| . Especially the system with = r (r = v/ ^~ 1 ) has been extensively studied 113,17 . 



The spectrum around E = shows a self-similar structure and the clear scaling behavior is 
observed. The wavefunctions are critical and some of them show the multifractal behavior 

In Ref. |0] , it is shown that the Hall conductivity carried is quantized to be an integer in 
units of e 2 /h when the Fermi energy Ep is in a gap. This integer is given by the first Chern 
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number of the fiber bundle on the magnetic Brillouin zone [f[0|j . It is the total vorticity of 



the U(l) phase of the Bloch wavef unctions which is a topological invariant. On the other 
hand, the conductivity of the system is also described by the edge states |y],[12|. It also 
has a topological origin. The winding number of the edge states in the complex energy 
plane gives the Hall conductivity |T2|. Dougot and Stamp discussed AC conductivities for 



a commensurate flux = p/q when a particle density p = <\> ||28|| . In this case, the Fermi 
energy is in the largest gap and the system is an insulator [p0j , [2ll1 . 

In this paper, we consider both longitudinal and transverse (Hall) conductivity. Espe- 
cially we are interested in the incommensurate flux limit. When the flux cf) is irrational, 
there are infinite number of bands with zero width (Cantor set). When the Fermi energy 
Ep is not in a gap with a finite width in the Cantor set, it is highly nontrivial whether the 
system is metallic or not. 

In Sec. ||, we derive the expression for the conductivity o ^^{uS) for a commensuarate flux 
<p = p/q when Ep is at an arbitrary position in the spectrum using the Kubo formula. In the 
following sections, the longitudinal conductivity is discussed in details. In Sec. JTTJ, systems 
with a sequence of rational numbers {<pi} are treated numerically in order to understand the 
system in the incommensurate flux limit. In Sec. [TV], we discuss the incommensurate limit 
by taking an appropriate scaling argument. Sec. is a summary. 

II. DERIVATION OF THE CONDUCTIVITY 

In this section, we derive an expression for the conductivity <j^ v ^{uj) for a rational flux 
= p/q by the Kubo formula [p9j] . The real part ^.cr xx ^{ijj) will be discussed in details. 
Let us rewrite the Hamiltonian (II . II) as 



^ = -E(^i^ + e-^ c t c ^). (2.2) 



where I denotes a two-dimensional lattice site and v is a unit vector along the v direction. 



The phase factor 9* is on the link between sites I and I + v. We consider a response of the 
system to a uniform time-dependent electric field E(t) = (E x (t) , E y (t)) . Since the electric 
field E%j on a link (ij) is given by E^ = — | -^j-, the Hamiltonian H with .E(t) is given by 
changing the phase factor 9\ in Eq. (|]|) to 9\ + f A v (t) with = Up to the 

first order in E(t), the Hamiltonian is expanded as 

H = H° + H'(t), H'{t) = -Y,J°vA»{t), (2.3) 

V 

where 

^ = '1^^-^'^, (2.4) 

I 

which represents an unperturbed current operator along the /x direction. The current oper- 
ator is 

f 

It is expanded as 

J„ = J° + ^H°A»(t), (2.6) 

up to the first order in E(t). By the Kubo formula, the induced electric current 5(J M (i)) is 
given by 

6(J»(t)) = ¥ (Hl)A»(t) + -Y: I dT([J°(T),J°})A»(t-r), (2.7) 
Ji(r) = e^jy^. (2.8) 

where (■ • •) denotes the ground state expectation value. By Fourier transformation, we get 

.2 n ! 



aio; — iO 



£», (2.9) 



with 



K lv {w)= / rfr([jJ(r),J°])e— . (2.10) 

J 



Let us take the Landau gauge, that is, 9~ = 0, 9~ = 27r0m where / = (m, n). In a momentum 
representation, we have 

r fl k ~ n 

H ' ' = L,,TM' cHk)H ' (kHk) - (2 - U) 



with 



H\k) = J2H l/ (k), (2.13) 

V 

c\k) = [c\k x + 2tt0, *g , ■ ■ • , c f ( fc x + 2tt0 g, fcj] , (2. 14) 

{H Jk)}ij = -2cos(^ + 27T0J) (J ij3 ., (2.15) 



.0 , 



{H y (k)} itS = -(e-^5 iJ+ i + e^^-x), (2.16) 

{««(*!) }.j = 2sin(/c c + 2vr0j) S id , (2.17) 

{^(fc)}^ = i^Hj+i - e'Hj-i), (1 <i,j< q), (2.18) 

where = (k x ,k y ) is defined in the magnetic Brillouin zone (M.B.Z.) (— ir/q < k x < 
n/q, —tc < k y < 7r), and the operator d{k x , k y ) is defined by c\k x , k y ) = Y, m ,n e l( * kxm+kyn) c^^. 
Let us discuss a one-particle state of the Z-th band 

\¥) = f £L jz^{k)c\k x + 2^k y m. (2.19) 

The Schrodinger equation H°\w) = E \ w) is reduced to 

H°(k)\k,l) = E l {k)\k,l), (2.20) 

with 

\k,l) = [Vl(fe), • ■ ■ , 1> l q (k)]\ (k, l'\k, l) = 5v, h (2.21) 
From Eqs. ( |2.12[) and (|2.10| ), we have 

K ( w ) = y^J_ / d 2 A;V^'(fc)^(fe)- /' (2.22) 

v^(k)^(k,l\vJk)\k,l'). (2.23) 



Similarly we have 



e l Jk) = (k,l\H°(k)\k,l), 



(2.24) 
(2.25) 



where V is the volume of the system and f l (k) = 6(Ep — E l (k)). Since 5(J fl (uj)) 
V T T, v v»v,<I>{u)E v {<jj), we have 



2ni jm.b.z 



d 2 kJ2e l ,(k)f(k) 



1 



hw-io ^ 

f(k)-f(k) 



l 



+ 2m Lb.z. ^ V " {k)vl ^ k) hu - iO + EV{k) - E\k) hu-iQ- 



(2.26) 



Let us rewrite the second term as 
1 



2m Jm.b.z. Y vti L E l (k) 



f(k) - f\k) f 1 
;*(fc)Uu;- 



1 



2tt/ ./m.b.z. ^' J £''(*;) -£7'(]fe) - iO hco-iO + E l '(k) - E l {k) 



^.„, E l '(k) E l '(k)-E l (k) 



1 



Tio; — zO 



27TC ./m.B.Z ^ " 



Using the formulae 



and 



<9 



(k,l\ w \k,l f ) 



E l {k) - E l '{k) Hu-iO + E l '{k) - E l {k) 
<(fc) 



<(fc) 



E l '(k)-E l (k) 
dE l (k) 



dk„ 



. (2.27) 



(2.28) 



(2.29) 



the first term of Eq. ( |2.27 ) is written 



2m ./m.b.z. Y vti 
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< (fe)(fe,i'l^-|fe,z> - (M^MKW 



1 



Tia; — iO 



— / rf 2 fc 

27T2 JM.B.Z. 



X 



£/'(*) 



i 

1 



d 







dk„ 



2m Jm.b.z. 

27T2 JM.B.Z. 



(k,l\ 



[ dk„ 
dvjje) 
dk.. 



(k, l\v u (k) — \k, 1} - (k, — ( v u (k)\k, 1} ) + T^{k 1 l\v v (k)\k, 1} 







k, l) + 



dv l J(k) 



dk„ 



dk,j 
1 



1 



hu — iO 



huo — iO 



hu — iO 2m' 



vtWm.b.z. Y 



dE ^hcu-iO' 



(2.30) 
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Next let us separate the cu-dependent part in the second term of Eq. ( |2.27| ) and we have 

— ( AV/^W'^^fc^ - 

2m Jm.b.z. f£, " K ' vK ' E l (k)-E l \k)huj-i$ + E l \k)-E\k) 



2ni Jm.b.z. frf, '' 



(E l (k) - E l '(k)) 2 
f{k)-f\k) 



(E l (k) - E l '(k)) 2 haj-iO + E l '(k) - E l (k) ' 



(2.31) 



Using Eq. (ggg) and (k,l\^\k,l') = -lJL(k,l\\\k,l'), the first term of Eq. (|2T3lD is 



written 



2m jm.b.z. 

1 
2tH 

1 



d 2 * £/'(*)£ 



^'(fc)^(fc)-<(fe)<(fc) 



I'^l 



(E l (k)- E l '(k)) 2 
d .,. „. 9 



27U JM.B.Z. 



d 2 k-£f(k) 



^-/ M . R / 2 *i:/'wvx( fc , ( |v| fc , ( )^ 



(2.32) 



From Eqs. (|2~26|) , (|2T30|) , (ggg ) and ( ^32|) with df(k)/dE = -8{E F - E l (k)), we have 

— / d 2 A;V \v l Uk)\ 2 5(E F - E\k))— — 8„ v 

2m Jm.b.z. Y JJ tiw-iO * 

d 2 k^2f l (k)Vx(k,l\V\k,l)e„ 
i 

f(k)-f(k) 



2m Jm.b.z. 
1 



+ 



27T2 Jm.b.z. 



^^(fc)^) 



{E l {k) - E l '(k)) 2 Hu-iO + E l \k) - E l {k) 



. (2.33) 



For the Hall conductivity a xy ^(u), we have 



I, ,\ _TKNN | S (, \ 



TKNN 



1 



(2.34) 
(2.35) 



«*» = li L, ^ g ^W<'W ( ig^[g )2fa _ i0 + ^ fc) _ g , (fcr (2-36) 
Since erf Au = 0) = 0, the static Hall conductivity is given by cr™^ N which is a topological 



invariant in the M.B.Z. 
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For the longitudinal conductivity cr xx ^(u), we have 
<WM = <rg t<l> (u;) + (2.37) 

^^L^p^i^-^^ (2 - 38) 

<W"J - ^ J m b z d K 2, I « x (fc) I _ ^ /(fc))2 ^ _ . + EV{k) _ El{k) , 

where crg. ^co) is the so-called Drude term and <rf ^(u;) is the contribution from interband 
scattering processes. Real parts of them are evaluated as 

= D(q)5(hw), (2.40) 
= \ f d 2 k^\vm\ 2 S(E F -E l (k)), 

2 JM.B.Z. j 

| §E™(k)=E F dk \ V x im ( k )\ 2 ^ E \k)\ i E F is in the ™" th band ) 

(2.41) 

(Ejr is in the energy gap), 



and 



= \ I I v l x l '(k) \ 2 {f(k) - f(k)}6(hu, + E l '(k) - E\k))±- 

M' i 



~ 2^^, E g F ^,i W -^'(^^ k " (fc)l ME*W-E"(k))\' (2 ' 42) 

The Drude term ^t<T^ x ^(u) is a delta function at uo = with the weight being determined 
by the average velocity of the electrons at the Fermi lines (there are q Fermi lines). The 
second term 9fJcrf ^(u;) comes from the interband scattering processes between states (fc, V) 
and (fe, /). 

For a rational flux, the onset of 3?crf x ^(a;) is given by the minimum energy gap. When 
the Fermi energy is in the energy band, the Drude term 3?0"£. ^(cj) can also contribute. 

III. NUMERICAL RESULTS FOR THE IRRATIONAL FLUX 



In order to understand the incommensurate flux case, we approximate the irrational (j) 
by a series of rational fluxes which tends to the irrational as q — > oo. For a large q, the onset 
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of the $ta xx Auj) becomes small and it is difficult to distinguish $tcr xx ^(uj) from $ta xx ±(uj). 
Thus careful considerations are necessary to understand the behavior of filer xx ^(uj). 

Let us consider the systems with (pi = Fi/Fi + \ where F\ is the Fibonacci number defined 
by Fi + i = Fi + with Fq = F\ = 1. The flux <pi tends to <p = (the inverse of the 

golden mean for / — > oo). Calculations were done for I = 5 ~ 12 which correspond to the 
fluxes 2T) §g) §§ 5 ^H' 5§3 anc ^ IP* ^ ne one " dimensional integrations in Eq. fl2.41|) and 
Eq. ([2.42 ) were numerically performed by the double exponential formula. 

First Ep was put in an energy band. We kept Ep near zero since there is a energy band 
around E = if q is odd. The cases with even q were avoided since the density of states 



disappear linearly at E = [p5| -|2"7j|. The Drude weights D(q) are shown in Fig.l. We treated 
three systems with Ep = 0, Ep = — 0.475A/ and Ep = — 0.4995A; where A; is a width of 
the central band in the l-th stage. In some cases, we have data for <fi = 610/987(7 = 14). 
Since it is known that the spectrum shows a self-similar structure for every three fluxes 
((pi, (fii+3, ■ ■ ■) |0|, we fitted the results within these groups. Numerical results are well fitted 
by the following scaling form 

D(qi)=D-^, a = 0.828 ±0.002. (3.1) 
<?f 

where the constant D depends on the position of Ep but a is universal (see Fig.l). We 
calculated the interband scattering term fitcr xx ± (u) when Ep = and they are shown in 
Figs. 2. In the qi — > oo limit, the energy gap near E = becomes infinitesimally small and 
the interband scattering term fitcr xx ± (uj) can also contribute in uj ~ 0. Next let us put Ep in 
the energy gap near E = 0. We have calculated fitcr xx ,<^(^) when Ep is in the gap just above 
E = 0. Numerical calculations for even q are also included. In this case, the Drude term 
does not contribute and we have calculated only $lcr xx j, (uj) . Numerical results are shown in 
Figs. 3. 

As shown in Figs. 2 and Figs. 3, the self-similar structure is reflected on the u-dependence 
of fitcr xx l p l (uj) . This is consistent with the the self-similar behavior observed in the energy 
spectrum. The supports of fila^ ^u) vanish when the flux approaches to the irrational 
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value. 



IV. DISCUSSIONS 

Since ^tUxx^X^) shows the complex behavior, we analyze the numerical results in Sec. 
|nj carefully and discuss scaling behaviors of the longitudinal conductivity. Since 0; can 
be considered as i-th approximation of the irrational flux <p = t, the energy resolution of 
the calculations should be determined by the spectrum in the Z-th approximation. Thus 
^• a xx'<t>, ( u ) i s averaged over the energy window [u - Wi/2,uj + Wi/2] to define fta^iu) 

where W\ is determined by the energy resolution of the energy spectrum. We put uj = 
to consider the small uj behavior. Denote 3?cr^'?(0) by 5Rcr^' 5 ^(0/). We take Wi the same 
as the energy scale AEi of the energy spectrum around E = 0. Let us take AEi to be the 
energy gap around E = 0. The scaling form of AEi is numerically given by 

coiisti 

A Ei = — =-^, 77 = 1.8285 ± 0.0006. (4.2) 



This scaling behavior was already discussed in Ref. [15]. We choose the window W\ as 



Wt = ^, (4.3) 
Qi 

where Wq is a q- independent constant. The number of bands M within the window Wi 
is determined by W which is independent of /. The average conductivities 9?cx ave (0j) were 
calculated for three values of Wq which give M = 3,7, and 11 respectively. We consider the 
two cases separately. One is when the Fermi energy is in the energy gap and the other is 
when the Fermi energy is in the energy band. First let us consider the former case. The 
averaged values of 9ft<7^ e (0z) are shown in Fig.4.(a). The results suggest strongly the scaling 
behavior 

MvlM = c{W ) qf s , 5s = 0.99 ± 0.02. (4.4) 
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The coefficient c(Wo) depends on the system but it is independent of qi. Next let us consider 
the latter case when E F is in the band. We calculated 9fcr^(0z) similarly. The results 
are shown in Fig.4.(b). They also suggest the scaling form Eq. ( [4.4| ) with the exponent 
5s = 0.99 ± 0.07. In this case, the Drude term $la xx aXoS) = D(q)5(hu) also contributes to 
9^0ave(0z) since the Fermi energy is in the energy band. The contribution from the Drude 
term is obtained as 

1 r w i/ 2 

^ZMi) = T7F / , 3*<^0 du' = c'(W ) q D , (4.5) 

Wi J-Wi/2 

where (/(Wo) is independent of q\. The scaling exponent a for the Drude weight is defined 
from Eqs. (JO) and (O) to be 



a = r] — 5]j. (4.6) 

The value of the exponent 5d is identical to 5 s within our numerical accuracy. 

Define 5 = 5s — 5^, then the average conductivity 9far ave (0/) has the scaling form 

3fo ave (0O = c(W ) qf . (4.7) 

This form is universal in the sense that the exponent 5 is independent of the position of 
Ep and the way the approximation of the irrational flux is made. Furthermore, it does not 
depend on whether the Fermi energy is in an energy gap or not. From Eqs. Q and (§77|), 
we have the scaling 

a^wM = for <*> ( 4 - 8 ) 

where 

7 = - ~ 0.55. (4.9) 
V 

The naive argument would lead to ^tcr XXjl p(uj) ~ c5(u) + d jus when Ep is in the energy 
band (metal) and ^to xx ^{uS) ~ c'/u when E F is in the energy gap (insulator) as u tends 
to zero. Our result Eq. (f4.8|) is different from both of these simple expectations. We shall 
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give the correct analysis below. Let us assume that the matrix elements in Eqs. ( |2.41|) and 
( |2.42j ) to be constant \v% (k)\ = v x . Then the Drude weight D(q) can be estimated as 



D{q)^\v x 2 i dk I =2^N m (E F )v x 2 . (4.10) 

where N m (E) = A § E m( kx ^ =E dk | V ^ m | is a density of states of the m-th band. Since 
there are q bands in the total spectrum, integrated density of states for each band is 1/q 
and N(Ep) ~ ~^E' wnere AE, the width of the band including E — 0, is of the order of 
AE ~ q~ n (see Eq. (|4.2| )). Also we may estimate v x = AE/Ak x where Ak x ~ 1/q is a size 
of the magnetic Brillouin zone in the x direction. Thus we have 

D(q) ~ const, x q 1 '" 1 . (4.11) 

It suggests a scaling relation 

a = rj - 1, (4.12) 

which is consistent with the numerical result of a = 0.828 ± 0.002 and 7] = 1.8285 ± 0.0006. 
It gives the scaling relation between the Drude weight and the spectrum. Note the scaling 
index for the Drude weight is a and the scaling index for the spectrum is r]. From Eq. ( |4.6|) 
it implies 



>D 



1. (4.13) 



Next let us consider the contribution from the interband scattering Eq. ( |2.42|) . Since the 
onset of ^ta^x 6, i s °f the order of AE, we have 

^L(^) ^ ^v x 2 N F N c , (4.14) 

where A^ c (~ 0(1)) is a number of a possible combination of the energy bands which can 
contribute the process. Thus we can estimate 3?<j^ e (0/) as K(rf ve (0i) ~ XEil^E) 2 -^^ = q. 
It suggests 

Ss = I- (4.15) 
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V. SUMMARY 



We have considered both of the longitudinal and transverse conductivities a^^cu) of 
the two-dimensional electrons on the square lattice at an arbitrary filling factor. Especially 
the longitudinal conductivity <7 xx ^{u) is discussed in details. It consists of two terms: the 
Drude term cr^ ^u) and the interband scattering term <rf ^(u;). 

We numerically calculated the real parts 3?cr£. ^(a;) and ^ta^^uo) for a series of sys- 
tems with rational fluxes {4>{\ which converge to the irrational value r = i^ 1 ■ The case 
when the Fermi energy E F lies near E = was treated in details. Two parts of the longi- 
tudinal conductivities Sftcr.^^ (a;) and 3?cr^ x ^ (a;) show the same scaling behavior in the in- 
commensurate limit. Taking into account the self-similar structure of the energy spectrum, 
we performed the averaging procedure near uj — 0. It revealed the clear scaling behavior 
$ta xx (u;) ~ 1/uo 1 ~ 0.55, uj — > 0). This result is quite distinct from the commensurate 
case. 
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FIGURES 

FIG. 1. Drude weights D (q) for systems with three Fermi energies Ep in the central band; (a) 
E F = 0, = f , if, §, and §f, (b) E F = 0, = f|, iff, and §f , (c) £ F = -0.475A,, = f , 
if, f, and §f, (d) E F = -0.475A,, = ± §f , iff, and §f, (e) £ F = -0.4995A,, = § , if , 
|| , and |ff , and (f) E F = -0.4995 A/, = ^, §§, |||, and §f where A/ is a width of the central 
band. 

FIG. 2. Kerf., (a;) for systems with rational fluxes {0/} with Ep = 0; (a) 5 = 8/13, (b) 

6 = 13/21, (c) 8 = 34/55, (d) 9 = 55/89, (e) 0n = 144/233 and (f) i2 = 233/377. Note the 
differences in the energy scales. 

FIG. 3. Kof^H for systems with rational fluxes {<&}; (a) 5 = 8/13, (b) 6 = 13/21, (c) 

07 = 21/34, (d) 8 = 34/55, (e) 9 = 55/89, (f) 0i O = 89/144, (g) 0n = 144/233 and (h) 
012 = 233/377. The Fermi energy Ep for each of them lies in the gap just above E = 0. Note the 
differences in the energy scales. 

FIG. 4. Averaged conductivities 3?cjf ve (0i) in two cases; (a) when Ep is in the energy gap, 
and (b) when Ep is in the energy band. ^<rf ve (0/) was calculated for three values of Wq which 
give M = 3,7, and 11, respectively; (a) (1)M = 3, = ^, §, iff, (2)M = 3, = if , §, |ff, 
(3)M = 3, = ||, ff, (4)M = 7, = §, iff, (5)M = 7, = if , §, |ff , (6)M = 7, = §, 
^, (7)M = 11, = ^, §, iff, (8)M = 11, = if, §§, |ff, and (9)M = 11, = §, ^, and 

(b) (1)M = 3, = A, 34 ; 144^ (2)M = 3> ^ = 13 j K 233^ (3)M = 7>0 = ^ j* 144 ? (4)M = 7> 

</> = i> §> ( 5 ) M = 11,^ = & 1> li' and (6)M = 11, = if, §, |ff . 
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